Electron-Phonon Coupling and the Soft Phonon Mode in TiSe2 
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We report high-resolution inelastic x-ray measurements of the soft phonon mode in the charge- 
density- wave compound TiSe2. We observe a complete softening of a transverse optic phonon at the 
L point, i.e. q = (0.5, 0, 0.5), at T « Tcdw- Renormalized phonon energies are observed over a large 
wavevector range (0.3, 0, 0.5) < q < (0.5, 0, 0.5). Detailed ab-initio calculations for the electronic and 
lattice dynamical properties of TiSe2 are in quantitative agreement with experimental frequencies 
for the phonon branch involving the soft mode. The observed broad range of renormalized phonon 
frequencies is directly related to a broad peak in the electronic susceptibility stabilizing the charge- 
density-wave ordered state. Our analysis demonstrates that a conventional electron-phonon coupling 
mechanism can explain a structural instability and the charge-density- wave order in TiSe2 although 
other mechanisms might further boost the transition temperature. 

PACS numbers: 71.45.Lr, 63.20.kd, 63.20.dd, 63.20.dk 



The origin of charge-density- wave {CDW) order, i.e., 
a periodic modulation of the electronic density, is a 
long-standing problem relevant to a number of impor- 
tant issues in condensed matter physics, such as the 
role of stripes in cuprates [l| and charge fluctuations in 
the colossal magnetoresistive manganites [2]. Chan and 
Heine derived the criterion for a stable CDW phase with 
a modulation wavevector q as Q 
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where rjq is the electron-phonon coupling (EPC) ma- 
trix element associated with a mode at an unrenormal- 
ized energy of ujtare, Xq is the dielectric response of the 
conduction electrons, and Uq and Vq are their Coulomb 
and exchange interactions. Static CDW order typically 
is taken as a result of a divergent electronic susceptibility 
Xq due to nesting, i.e. parallel sheets of the Fermi sur- 
face (FS) separated by twice the Fermi wavevector 2k f. 
Electron-phonon coupling (EPC) is required to stabilize 
the structural distortion and, hence, an acoustic phonon 
mode at the CDW wavevector qcDW = 2A;/ softens to 
zero energy at the transition temperature Tc dw 0; 0! ■ 
However, when electronic probes reported only small and 
not well nested Fermi surfaces this scenario has been dis- 
carded for the prototypical CDW compound TiSe2 0-01 ■ 

Alternative scenarios such as indirect or band-type 
Jahn- Teller effects |6|-|8| and, most prominently, exci- 
ton formation [1, 0, lid ] are discussed in the theoretical 
as well as experimental literature. More recently, van 
Wezel et al. have invoked a model including exciton 
formation as well as EPC JJL] and have shown that it 
can explain data from angle- resolved photoemission spec- 
troscopy (ARPES) formerly taken as evidence of an 



excitonic insulating phase in TiSe2 0^- Determining the 
origin of CDW formation in TiSe2 is all the more im- 
portant with respect to the nature of superconductivity, 
which emerges both as function of pressure [l^l and Cu 
intercalation [3]. In particular, pressure induced super- 
conductivity is expected to be closely linked to the nature 
of the parent CDW state. 

In this letter, we report a temperature dependent 
study of the CDW soft phonon mode in TiSe2 by high- 
resolution inelastic x-ray (IXS) scattering. We inves- 
tigated both phonon energies as well as the phonon 
linewidths. We observe a complete softening of a trans- 
verse phonon branch at the L point close to the CDW 
transition temperature Tcdw- Calculations of the lat- 
tice dynamical properties based on density-functional- 
perturbation theory (DFPT) describe our results very 
well and demonstrate that conventional EPC can stabi- 
lize CDW order contrary to the current discussion in the 
literature 0,i,[lii. 

Our sample was a high-quality single crystal with di- 
mensions (3 X 3 X 0.05) mm'^ that shows a sharp onset 
of a (2 X 2 X 2) CDW superstructure at Tcdw ~ 200 K 
as observed by synchrotron x-ray diffraction (XRD) fl^ . 
The high-resolution IXS experiments were carried out at 
the XOR 3-ID HERIX beamline of the Advanced Pho- 
ton Source, Argonne National Laboratory. The inci- 
dent energy was 21.657 keV and the horizontally scat- 
tered beam was analyzed by a spherically curved sili- 
con analyser (Reflection [18,6,0]). The full width at 
half maximum (FWHM) of the energy and wavevec- 
tor space resolution was 2.4 meV and 0.07 A~^, respec- 
tively. The components {Qh,Qk,Qi) of the scattering 
vector are expressed in reciprocal lattice units (r.l.u.) 
{Qh,Qk,Qi) = {h * 27r/a, fc * 27T/a,l * 27r/c); with the 
lattice constants a = 3.54 A and c = 6.01 A of the hexag- 
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FIG. 1. (Color online) Temperature dependence of soft 
phonon mode at Q = (—0.5, 1,0.5) {L point), (a)-(c) Energy 
scans for temperatures 204 K < T < 300 K. Solid (red) lines 
are fits consisting of a damped harmonic oscillator (inelastic) 
and a pseudo-voigt function (elastic) (blue dashed lines), (d) 
Soft mode energy ujq observed via IXS as function of temper- 
ature. The dashed line in (f) is the temperature dependence 
deduced from thermal diffuse scattering reproduced from Ref. 
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FIG. 2. (Color onhne) Energy scans at Q = (— /i, 1,0.5), h = 
0.35 - 0.5, for temperature T = 190K (« Tcdw). Solid lines 
represent the total fit result consisting of a damped harmonic 
oscillator for the inelastic and a pseudo-voigt function for the 
elastic scattering (blue dashed lines). The inset of (d) shows 
the scan &t h = 0.5, which can be well described just by the 
pseudo-voigt function as used in (a)-(d). 



onal unit cell (PSml, 164). 

Measured energy spectra were fitted using a pseudo- 
Voigt function for the elastic line with a variable ampli- 
tude and fixed lineshape determined from scans through 
the CDW superlattice peak at base temperature. The 
phonon was fitted by a damped harmonic oscillator 
(DHO) function, where the energy uiq of the damped 

phonon is given by LOq — Juj^— [l^l, where ujq is the 



phonon energy renormahzed only by the real part of the 
susceptibility, and F denotes the phonon linewidth, which 
is closely related to the imaginary part of the suscepti- 
bility. The situation, where the phonon becomes over- 
damped and freezes into a periodic lattice distortion, is 
defined when the damping ratio F/cjg reaches its critical 
value, F/wq = 1. 

Density-functional-theory calculations were performed 
in the framework of the mixed basis pseudopotential 
method [l^. The exchange-correlation functional was 
treated in the local-density approximation (LDA) . Norm- 
conserving pseudo-potentials for Ti and Se were con- 
structed with Ti 3s and 3p semicore states included in the 
valence space. Phonon frequencies and electronic contri- 
butions to the phonon linewidth were calculated using 
the linear response technique or density functional per- 
turbation theory (DFPT) [3l in combination with the 
mixed-basis pseudopotential method ^] . To resolve fine 
features related to the Fermi surface geometry, Brillouin- 
zone (BZ) integrations were performed with a dense 
hexagonanal 24 x 24 x 8 fc-point mesh (244 points in the ir- 



reducible BZ). The standard smearing technique was em- 
ployed with a Gaussian broadening of cr = 0.15 eV. Dis- 
played results were obtained for the undistorted hexag- 
onal structure using experimental lattice constants and 
optimized internal parameters. The calculated electronic 
structure and lattice dynamical properties are in good 
agreement with previous results from LDA 2^, 21 [ . 



Fig. [T] shows the temperature dependence of the soft 
phonon mode measured at the i^point, i.e. Q — 
(—0.5,1,0.5), from temperatures just above Tcdw up 
to room temperature (Fig. [T^-c). At room tempera- 
ture we find a phonon energy ix)q — 5.3 meV in good 
agreement with previous neutron scattering data |22| . 
On lowering the temperature, softens to zero energy 
(within the experimental error) already at T = 204 K. 
At Tcdw ^T < 204 K, the phonon is overdamped and 
can only be observed as a broad quasi-elastic peak un- 
derneath the sharp elastic line, which becomes gradually 
more dominant below 200 K (not shown). For T < 190 K 
low energy phonon scattering directly at the L point 
was undetectable due to the strong CDW superstructure 
peak. Overall, our results roughly agree with energies 
of the soft mode deduced from thermal diffuse scattering 

We studied in detail the dispersion of the transverse 
optic (TO) phonon starting at 17meV at A and ending 
in the CDW related soft mode at L, i.e. q = (— /i, 0, 0.5) 
and < /i < 0.475 [? ]. We note that at h ^ 0.33 there 
is an anti-crossing of this TO branch with an acoustic- 
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FIG. 3. (Color online) (a) Dispersion of the soft phonon mode 
at r = 190 K (f» TcDW, filled symbols) and T = 30K (open 
circles) along Q — {—h, 1, 0.5), h — 0—0.475. Dashed lines are 
results of DFPT calculations using the experimental crystal 
structure and a — 0.15 eV. The solid (blue) line indicates the 
dispersion of the soft mode character across the anti-crossing 
at h St 0.33. Frequencies of the upward dispersing branch 
(e.g. see fig. 2(a)) agree well with the calculation but are not 
shown for the sake of clarity, (b) Damping ratio of the DHO 
at T = 190 K (filled symbols) and T = 30K (open circles). 



like upward dispersing branch. The character of the soft 
mode is transferred to the energetically lower branch for 
h > 0.33. For simplicity, we speak of the soft-mode 
dispersion as represented by the dashed (blue) line in 
Fig. [3^, which involves parts of the TO branch {h < 0.33) 
and the acoustic-like branch {h > 0.33). A set of energy 
scans taken at T = 190 K is shown in Fig. [51 

At h = 0.35 (Fig. [2^), we observe two phonon peaks, 
where the one at lower energies corresponds to the char- 
acter of the soft mode starting at 17meV at A. The 
second peak near ITmeV represents the continuation of 
the upward dispersing acoustic-like branch in good agree- 
ment with structure factor calculations by DFPT. This 
mode does not show any significant temperature depen- 
dence. Results of its dispersion are not shown in Fig. [3] 
for the sake of clarity and will be discussed in detail else- 
where. 

The energy of the phonon peak related to the soft 
mode already decreases for h > 0.25 and, after the 
anti-crossing, shows a quasi-linear dispersion relation for 
h > 0.4 (Fig. [2j)-d) , which extrapolates to a zero phonon 
energy at ft. = 0.5 (Fig. |3^). At h = 0.5 however, the 
strong CDW superlattice reflection made it impossible 



to directly measure the quasielastic scattering expected 
for the critically damped phonon. The calculated disper- 
sions agree very well with these results. Corresponding 
to the strong renormalization of the phonon energy, we 
observe an increase of the damping ratio F/wg, which 
stays just below 1 for h = 0.475. But just like the en- 
ergy, the damping ratio also shows linear behavior for 
h > 0.4 and extrapolates to an overdamped phonon at 
h = 0.5 (Fig. [3b). 

The range in q over which phonon energies are renor- 
malized as function of temperature allows us to infer 
how sharp corresponding electronic signatures, e.g. in 
Xq in a nesting scenario, might be |J|. We note that, 
in principle, the q dependence of the EPC matrix ele- 
ment rjg can determine the periodicity of the low tem- 
perature superstructure as well [2j]. Because the soft 
mode at T = 300 K is still far from being unrenormalized 
{T/ujq{T = 300 K) = 0.5), we compare the dispersion at 
T = 190 K with results for T = 30 K in the CDW ordered 
state, as the latter give a good idea of the range of renor- 
malized phonons. Fig. [3| clearly shows, that the phonon 
energies over the range 0.3 < h < 0.5 are renormalized 
across the CDW transition. This large range in h is in 
contrast to compounds, where a sharp FS nesting drives 
the CDW formation and renormalized phonon energies 
are observed only over a much smaller q range 2^, 2q . 

Comparing the experimental results in more detail to 
the calculations, it is important to note that DFPT is 
a single particle response theory and, therefore, cannot 
include, e.g., exciton formation. However, EPC is cal- 
culated in great detail and, in particular, the electronic 
contribution to the phonon linewidth 2"f can be deduced 
(Fig. [3]). 27 of a specific phonon mode in turn can be 
approximately expressed as [27] 



27 = ATTLUbare X \r]g\'^ X J DOS 



(2) 



where the electronic joint density-of-states (JDOS) in 
turn is closely related to the imaginary part of the 
electronic susceptibility Xg- We see that the q range 
over which 27 shows a strong enhancement matches 
closely the range of observed renormalized phonon en- 
ergies (Fig. [3^). Furthermore, 27 also qualitatively ex- 
plains the q dependence of the observed linewidth F 
very well (see Fig. [Hd). Quantitatively, the calculated 
linewidth is a factor of 2.5 smaller than observed in ex- 
periment. Here, we note that similar observations have 
been made in the CDW compound NbSe? [3] and the 
conventional superconductor YNi2B2C [28|. In the lat- 
ter, it is evident that only the lowest acoustic modes, 
which exhibit strong EPC, show this effect of too small 
calculated linewidths, whereas the linewidths of higher 
energy modes are predicted in much better agreement 
[2^ . Therefore, we think that soft modes as in TiSe2 are 
subject to an increased anharmonicity close to the struc- 
tural phase transition refiected in an increased linewidth, 
which is not captured in our theory. 
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FIG. 4. (Color online) (a) Results of the calculation are shown 
as function of h along the A — L line, i.e. q = {h, 0, 0.5) with 
h = — 0.5: (1) the calculated electronic contribution to the 
phonon linewidth 2j of the soft phonon mode (FWHM, solid 
line, left hand scale), (2) the total electronic JDOS (open 
symbols, right-hand scale) and (3) interband contributions to 
the JDOS between the Se 4p and Ti 3d bands near V and L 
(filled symbols, right-hand scale), (h) Experimental linewidth 
r (dots) compared to scaled 27 (x 2. 5-1- off set, solid line). 



We now discuss the different contributions to 27, i.e. 
\rjq\ and the JDOS, according to equ. [H The fact that 
27 for the soft phonon does not follow the increase of the 
total JDOS below h = 0.2 is due to the presence of dif- 
ferent bands near the Fermi energy. Whereas the JDOS 
at small h is dominated by scattering processes within 
the Se Ap bands around the F point, at large h scattering 
between Se Ap and Ti id bands dominate (Fig.|3ti). Ap- 
parently, rjq of the soft mode only samples the latter scat- 
tering processes. Using equation [21 ujbare — 17meV and 
the partial JDOS of Se Ap - Ti Sd scattering processes, 
we can calculate \'qq\. However, this is only possible for 
h > 0.25 as scattering processes of the same type seem 
to contribute to the EPC of different phonon modes for 
h < 0.25 and a clear determination of the relevant JDOS 
is not possible anymore. We calculate \T]q\ = 0.098 meV 
a.t h = 0.5 and a 33% decrease towards h = 0.25. For the 
same q range, the partial JDOS of Se Ap - Ti 3d scatter- 
ing processes decreases by a factor of six. Hence, the q 
dependence of the JDOS is the dominant factor in the 
large increase of the phonon linewidth. The correspond- 
ing phonon renormalization then leads to the structural 
instability at the L point. 

Our results explain the CDW phase transition in TiSe2 
in the conventional picture of an increased electronic sus- 



ceptibility at q_cDW aided by a strong EPC matrix ele- 
ment. This is not in contradiction to ARPES data, on 
the basis of which FS nesting previously has been ruled 
out as a driving force of the CDW transition in TiSe2 . 
Indeed, our results confirm that TiSe2 does not exhibit 
a divergent JDOS. However, the peak in the JDOS is 
connected with an enhancement of Re x 2l| , which can, 
according to equ. [U tip the balance between undistorted 
and CDW ordered phases if stabilized by a sufficiently 
large rjq. We note that DFPT is a zero temperature 
approach and, therefore, cannot calculate a transition 
temperature corresponding to, e.g., the calculated EPC 
strength. So, EPC might be only one ingredient for driv- 
ing the CDW transition in TiSe2 , but our analysis shows 
that its importance has been greatly underestimated. 

In conclusion we have reported a detailed analysis of 
the soft phonon behaviour in the CDW compound TiSe2 . 
The experimental observations are well reproduced by 
ab-initio calculations, which demonstrate that electron- 
phonon-coupling is strong enough to stabilize a structural 
phase transition in TiSe2 at finite temperatures. Our 
work emphasizes the essential part of electron-phonon 
coupling matrix elements and provides a new route to 
seemingly mysterious phase transitions as far as the elec- 
tronic footprint is concerned. 
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